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1. Introduction 

Quark-Gluon-Plasma (QGP) is a state of matter in the Strong Interaction system which is 
actively sought through theoretical study and heavy ion collision at high energy (or high temper- 
ature) and at low energy (or high baryon density) and must have played an important role in the 
evolution of early universe [|I|]. Unfortunately at physical quark masses, chiral symmetry restora- 
tion/deconfinement transition in QCD from normal hadronic matter to QGP matter is a cross-over 
transition and lacks characteristic behaviors of order parameter in a first order or a second order 
thermodynamic phase transition across the thermodynamic transition. Thus, finding and under- 
standing physical phenomena which can serve as signatures for QGP matter formation is quite 
crucial. 

Physics of heavy quarkonium in medium may be one of such candidate phenomena In 
contrast to large uncertainty in our understanding of light hadron production and decays, inclu- 
sive productions and decays of heavy quarkonium in proton-proton or proton-anti-proton colli- 
sion, which are required for the comparison with those in heavy ion collisions, is well understood: 
Non-Relativistic QCD (NRQCD) factorization in which the scale above heavy quark (like charm 
quark or bottom quark) mass is "integrated away" and the smallness of heavy quark velocity at the 
quarkonium rest frame offers "small parameter" is applicable in heavy quarkonium Short dis- 
tance effects can be calculated via perturbative expansion in and long distance non-perturbative 
matrix elements can be computed either using lattice NRQCD at zero temperature [Q, ||, |6|] or can 
be fitted using experimental data (see e.g., [^). Based on our understanding of quarkonium physics 
in p-p collision, medium effect for quarkonium in heavy ion collision may be better disentangled. 
In particular, bottomonium system offers a better probe for QGP formation since there are clearer 
seperations between the long distance physics and the short distance physics, compared to those in 
charmonium system due to smaller quark velocity (v^ ~ 0.1 compared to ~ 0.3 for charmonium 
system) at the quarkonium rest frame. 

In recent years, we studied the behavior of S-wave and P-wave bottomonium states for a 
temperature range QATc <T <2.lTc using NRQCD in non-zero temperature on anisotropic lattices 
(12^ X Nt) with two light quark flavors [||, ^. Details on how dynamical gauge field is generated 
and what the simulation parameters are can be found in ^ |l^. Our calculation method is 
explained in [||, ^]. We expect that in-medium effect such as temperature (T) effect or baryon 
chemical potential (pi) effect can be accommodated by NRQCD as long as ^ <C 1 or ^ ^ 1 where 



M is the heavy quark mass (see, e.g. [11] for baryon chemical potential efffect in NRQC2D). 
Indeed our study [^, ^ meets such expectation. We found that P-wave bottomonium states (Xbj) 
melt immediately above Tc, and bound state signals for S-wave channel (t]&,T) survive above 
[Q]. Closer inspection of the spectral functions for S-wave correlators (;^/,,T), which is obtained 
by the Maximum Entropy Method (MEM) with NRQCD kernel, shows that (IS) peaks appears to 
survive and the excited state peaks are suppressed [^. This observation is consistent with recent 
CMS experiment which shows sequential suppression of (35), (2S) T states and survival of (1S)T 



[12, 13] 



Here, we discuss further lattice NRQCD study of ours on bottomonium at non-zero tempera- 



ture. First, bottomonium with non-zero momentum is studied (for full report, see [14|). Quarko 



nium production in hadron collision proceeds: parton evolution from hadron, then partonic scat- 
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Figure 1: Quarkonium production processes: (a)fragmentation from gluon (top), (b) recombination of ther- 
malized heavy quark-heavy anti-quark (bottom) 

tering process (gluon fusion, heavy quark recombination, gluon Compton scattering, and light 
quark-light quark scattering). Among these partonic processes, quarkonium production is mostly 
dominated by gluon fusion. Gluon subsequently fragments into quarkonium. In this case, frag- 
menting gluon will have large energy and momentum. If quarkonium is produced in thermal bath, 
then it will be moving relative to thermal bath. In contrast, if quarkonium is produced through 
recombination of thermalized heavy quark and heavy anti-quark (albeit less probable), produced 
quarkonium is at rest in thermal bath. Schematic view on these two processes is shown in Fig. 
[I[ Since a thermal bath defines a preferred rest frame, "moving" effect will induce ^Mv^ shift 
in the mass compared to the mass of quarkonium produced at rest. This shift may be larger than 
experimental resolution for muon detection and may be observable. 

Secondly, we start to investigate systematic eiTors associated with our computation in [Q, 0. 
As a first step, we investigate how 10 ~ 20% change in bottom quark mass affects previous obser- 
vations made in [||, ^ by calculating bottomonium correlators with Mas = 3.7,4.14 and 5.4 and 
looking at the coiTclator ratios. 

2. S-wave bottomonium moving in thermal bath 

We consider r\}, and T which are moving with regard to the thermal bath, i.e., S-wave bot- 
tomonium states with momentum. Non-zero momentum quarkonium state is obtained by injecting 
momentum to the source for NRQCD quark propagator computation and projecting out the in- 
jected momentum component from the quarkonium propagator which is obtained by combining 
the zero-momentum NRQCD quark propagator with the momentum injected NRQCD quark prop- 
agator. Momenta,^= ^(«i,n2,«3) where H= (0,0,0), (1,0,0), (1,1,0), (1,1,1), (2, 0,0), (2, 1,0), 
and (2,2,0) with L = 12, have been considered. Since Ma^ = 4.5 is used for the bottom quark 
mass in NRQCD propagator computation and Mj,(i5) = 9.391 GeV and Mx(i5) = 9.460 GeV [|l|], 
these quarkonium momentum range amount to quarkonium velocity, 0.07 ~ v ~ 0.18 with 
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Figure 2: T spectral functions on 12 x 24 for various momenta 
^2^2 _ 4^.^j 2 sin^ |-. From the correlators, spectral functions defined as 

G{T,p)= J^K{r,(o)p{co,p), K{T,(o)=e-''' (2.1) 

are computed using MEM. Fig. ^ shows a typical spectral function as a function of \p\. For each 
momentum, the spectral function exhibits a distinct shape. From this, the peak and width can be 
extracted and Fig. |3| show extracted peak values and width values which corresponds the T(15) 
state. To a leading order in p^ and T, we see that the 15 peak roughly follows 

AE{p,T) ~ AE{p = 0, r = 0) + ^ + <«s)T (2.2) 
In general, continuum NRQCD dispersion is 

E=Mo + — ^ + --- (2.3) 

2M2 SMl 

Since we use ^(v^) NRQCD lagrangian to compute bottomonium propagator , Mq has (^(v"*) error, 
M2 has ^(v^) error, and M4 ^(v°) respectively. Using the form of Mq + 4^, (sin(A:,/2)^)/2M2 - 
(4I,(sin()t;/2)2)2)/8M| to fit Nt = SOT correlator, we get 11.05 GeV for the kinetic mass (M2) 
of T(1S) state with Ma, = 4.5. The fitting form Mq + 4i:,(sin(A:,/2)2)/2M2 gives 11.12 GeV for 



M2. Thus Eq. (2.2) is suitable to describe our lattice result. On the other hand, the 15 width is 
independent of p^ and is just proportional to T. 

This lattice NRQCD finding can be compared with effective theory prediction. NRQCD bound 
state at rest in thermal bath below the melting temperature has been considered in [ |l6| ] and non- 
relativistic bound state which is moving with regai^d to thermal bath has been considered in QED 
(in the rest frame of the bound state, i.e., bound state in moving thermal bath) [[l7|]. Finding in ^ 
for the width can be summarized as 

r„ Vi-v^ ^_ fi+v\ ^ 1 2 



r. 2v log^TT^j-i-r + ^K) (2.4) 
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Figure 3: T(15) energy vs. various momenta from the first peak of the spectral functions at each T /Tc — 
2.09, 1.86, 1.68, 1.40, 1.20, 1.05, and 0.42 for M = 16,18,20,24,28,32, 80 respectively (left). Y(15) width 
vs. various momenta from the width of the first peak in spectral functions (right) 

where Fq is the width of the bound state not moving with regard to thermal bath and Fy is the width 
of the bound state moving with the velocity v in thermal bath. With < 0.04, "moving effect" in 
the width is beyond the precision level of our computation, which explains negligible momentum 
dependence of the width in the right figure of Fig. |3[ 

3. Mass dependence of bottomonium correlators 

Our previous bottomonium study was done with bottom quark mass, Ma^ = 4.5. The same 
NRQCD correlator calculation is repeated with Ma, = 3.7(-20%),4.14(-10%) and 5.4(+20%) 
in order to investigate how small change in bottom quark mass affects bottomonium spectrum and 
how important this effect is compared to non-zero temperature effect. Since lattice NRQCD is 
an effective theory and is consistent when Ma^ ~ 1, all these choices (Ma^- = 0.62,0.75,0.9 for 
Mas = 3,7,4.5,5.4 respectively) are expected to be valid. 

Lattice NRQCD quarkonium correlator behaves as Ae^'^^ when a bound state is formed and 
fitting T correlator calculated on A', = 80 configurations gives AEy{M, T ~ 0). Using experimental 
value for Mj, we determine the unknown constant for NRQCD spectrum by 

M^'^P = Eq{M) + AEy{M, T = 0) (3.1) 

If non-zero temperature quarkonium states has a pole like T(1S) state 

Mr{T)=Eo + AEr{M,T) =M-^'''^ + AEr{M,T)-AEr{M,T = 0). (3.2) 

So, the ratio of the non-zero temperature quarkonium correlator to the Nt = 80 correlator can reveal 
the heavy quark mass dependence AE{M,T) although the true ground state behaviors for Nt = 80 
coiTclators and for A/^t = 16 or 20 correlators will not set in at large T. Fig. ^ shows such ratios for 
T = 45S{Nt = 16) (MeV) and T = 408(A^, = 18) (MeV). 



5 



lattice NRQCD at non-zero T 



Seyong Kim 




Figure 4: T correlator ratio of G{x;N,) to G{z;N, = 80) for A^, = 16 (filled symbol) and 18 (open symbol) 

In Fig. 0, the comparison between the filled symbol and the empty symbol shows the tem- 
perature dependence for the same M, and the comparison between the different symbols shows 
M— dependence for the same T. The correlator ratio change due to 10% change in the two different 
T is larger than due to 10 ~ 20% change in the 4 different M. Since the small change in heavy 
quark mass is less important than that due to the temperature change, lattice NRQCD is valid as an 
effective theory and systematic error due to heavy quark mass mis-tuning is small. 

4. Conclusion 

Using NRQCD formaUsm in non-zero temperature, we computed momentum dependent (T]fe,T) 
correlators and calculated spectral functions using MEM via NRQCD kernel. From the peaks and 
the width of momentum dependent spectral function, we find that there are observable effect on 
the energy of S-wave state moving in thermal bath but no observable effect on the width of S-wave 
state moving in thermal bath for v^p^jj^^ ~ 0.04. 

Among various sources of systematic eiTors involved in our result, bottomonium correlators 
calculation is repeated with 20% variation in bottom quark mass. By studying correlator ratios. 
We find that the temperature effect is more important than the heavy quark mass effect in S-wave 
bottomonium at the temperature around a few T^, which suggests lattice NRQCD as an effective 
theory for bottomonium in non-zero temperature is a consistent theory. 

Improved tuning for bottom quark mass and detecting thermal deviations from the standard 
dispersion relation is in progress for Nf = 2+1 flavor system with smaller spatial lattice spacing 
and a larger spatial extent. 
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